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Abstract. High energy 7-rays from individual giant molecular clouds contain 
unique information about the hidden sites of acceleration of galactic cosmic rays, 
and provide a feasible method for study of propagation of cosmic rays in the Galactic 
Disk on scales < 100 pc. I discuss the spectral features of 7r°-decay 7-radiation from 
clouds/targets located in proximity of relatively young proton accelerators, and spec- 
ulate that such "accelerator +target" systems in our Galaxy can be responsible for 
a subset of unidentified EGRET sources. Also, I argue that the recent observations 
of high energy 7-rays from the Orion complex contain evidence that the level of the 
£f} • "sea" of galactic cosmic rays may differ significantly from the flux and the spectrum 

of local (directly detected) particles. 
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Q\ . 1. Introduction 
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The solution of the long-standing problem of the origin of galactic 
cosmic rays (GCRs) depends, to a large extent, on advances of 7-ray 
astronomy. The basic idea of this belief is straightforward and concerns 



o 

1-^ ■ both the acceleration and propagation aspects of the problem - while 



Oh 



the localized sources of 7-rays exhibit the sites of particle accelera- 



tion, the angular and spectral distributions of diffuse 7-radiation give 
■ information about the propagation of CRs in the Galactic Disk. 

The studies of the diffuse galactic 7-radiation by EGRET aboard 
the Compton GRO have already made a significant contribution to 
the knowledge of spatial distribution of low energy CRs in our Galaxy 
(Hunter et al., 1997). Furthermore, many famous representatives of 
different classes of potential particle accelerators like pulsars, shell- 
type supernova remnants, giant molecular clouds and associated with 
them star formation regions, are identified by EGRET as sources of > 
100 MeV radiation (Hartman et al., 1999). It should be noted, however, 
that most of the EGRET sources at low- and mid- galactic latitudes 
still do no have clear counterparts at other wavelengths (Gehrels et 
al., 2000). The next generation major space-based 7-ray detector, the 
GLAST (Gehrels and Michelson, 1999) with its superior flux sensitivity 
(see Fig. 1) and good angular resolution will be able to reveal the nature 
of these hot spots. GLAST offers a broad - up to 100 GeV for bright and 
flat-spectrum sources - spectral coverage, which ensures the great role of 
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Figure 1. The energy flux sensitivities of GLAST and future IACT arrays. For 
comparison, the broad-band 7-ray spectrum of the Crab Nebula is also shown. 

GeV is the domain of ground-based gamma-ray astronomy, and could 
be best explored by so-called imaging atmospheric Cherenkov telescope 
(IACT) systems (see e.g. Aharonian and Akerlof, 1997). The very high 
energy (VHE) region of 7-rays from 100 GeV to 100 TeV presents a 
special interest, because these 7-rays carry direct information about 
the sources responsible for formation of the most energetic part of the 
spectrum of GCRs extending to the so-called 'knee' around 10 15 eV. 
The range of flux sensitivities that could be achieved by future IACT 
systems is shown in Fig. 1 (Aharonian, 1999). It is seen that GLAST 
and future IACT arrays could probe point sources in a broad energy 
region from 0.1 GeV to 10 TeV at the energy flux level between 10~ 13 
and 10 -12 erg/cm 2 s. For comparison, the flux of the Crab Nebula, the 
best studied galactic 7-ray source, exceeds 10 -11 erg/cm 2 s both at TeV 
and GeV energies (see Fig. 1). Thus, all extraterrestrial sources of high 
energy 7-rays with fluxes as low as 10 mCrab should be detected by 
GLAST and/or future IACT arrays. In particular, these instruments 
should be able to reveal the sites of most powerful accelerators of CR 
protons in our Galaxy, more specifically to detect 7r°-decay 7-rays from 
shell-type SNRs (Drury et al., 2001) as well as from other possible 
non-thermal galactic source populations like pulsars, X-ray binaries, 
microquasars, etc. (see e.g. Aharonian, 1999). 

In this paper I do not intend to review the high energy 7-ray sources, 
but rather I shall focus on a specific class of objects - giant molecular 
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clouds (GMCs) located in the vicinity of CR accelerators. The existence 
of the CR accelerator by itself is a necessary but not sufficient condition 
for effective 7-ray production; obviously we need the second component, 
a target. The GMCs seem to be ideal objects to play that role in our 
Galaxy. These objects are connected with star formation regions which 
are believed to be the most probable (with or without SNRs) birth 
places of GCRs (Montmerle, 1979; Casse and Paul, 1980). 

2. Giant Molecular Clouds as tracers of cosmic rays 

One of the principal parameters which determine the 7-ray visibility of 
GMCs is M 5 /dl pc , where M 5 = M cl /1O 5 M is the diffuse mass of the 
cloud located at a distance c?k pc = d/1 kpc from the observer. Assuming 
that 7-rays are produced by interaction of CRs with the ambient gas, 



where q-2s{> E^) = g 7 (> -E 7 )/10 -25 (H - atom) _1 s _1 is the 7-ray 
emissivity. In a "passive" GMC, i.e. in a cloud submerged in the "sea" 
of GCRs, assuming that the level of the "sea" of GCRs is the same as 
the proton flux measured at the Earth, 



the 7-ray emissivity above 100 MeV is equal to (7-25(> 100 MeV) = 
1.53 77 (H— atom)~ 1 s~ 1 , where the parameter 77 ~ 1.5 takes into account 
the contribution of nuclei both in CRs and in the interstellar medium 
(Dermer, 1986). A "passive" cloud can be detected at the EGRET sensi- 
tivity level, if M^/d^, > 10, taking into account large angular size (typ- 
ically, several degree or more ) of objects with M 5 /d£ pc ~ 10. At energies 
E » lGeV, J(> E) ~ 1.5 x 10~ 13 (£7/1 TeV)- 1 - 75 (M 5 / &l pc ) ph/cm 2 s 
(Aharonian, 1991), thus even for M 5 /d^ pc ~ 10 (there are only several 
clouds in the Galaxy with such a large value) , detection of 7-rays from 
"passive" clouds above 10 GeV is an extremely difficult task both for 
GLAST and ground-based 7-ray detectors. 

Nevertheless, searches for very high energy 7-rays from GMCs present 
a certain interest because of possible existence of high-energy CR accel- 
erators nearby or inside GMCs (Issa and Wolfendale, 1981; Ormes et al., 
1988; Aharonian, 1995; Aharonian and Atoyan, 1996). Indeed, assuming 
that CRs with total energy W p injected by a source into the ISM during 
their propagation reach at the instant t a radius R(t), the mean energy 
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Figure 2. Temporal and spectral evolution of CR fluxes at different (fO pc, 30 pc, 
and fOO pc) distances from an impulsive proton accelerator. The power-law proton 
spectrum with a = 2.2 and total energy W p = 1 50 erg is supposed. The curves 1, 
2, 3, and 4 correspond to the age of the source t = 1 3 yr, 10 4 yr, 10 5 yr, and 10 6 yr , 
respectively. The energy-dependent diffusion coefficient D(E) with power- law index 
S — 0.5 is supposed. Figure 2a (on the left) and Figure 2b (on the right) are for 
Dio = f 28 and D w = f 26 cm /s, respectively. The hatched curve shows the local 
(directly measured) flux of CR protons. (From Aharonian and Atoyan, f996). 

density of CRs is equal to w p » 0.5(W P /10 50 erg)(R/100pc)~ 3 eV/cm 3 . 
Thus, in the regions up to 100 pc around the CR accelerators with 
W p ~ 10 50 erg, the density of relativistic particles at some stages, 
depending on the time history of particle injection and the character 
of their propagation in the ISM, may significantly exceed the average 
level of the "sea" of GCRs, w G CR ~ leV/cm 3 . In Fig. 2 the differ- 
ential fluxes of protons, at distances R =10, 30, and 100 pc from an 
"impulsive" accelerator, with total energy W p = 10 50 erg released in 
the form of relativistic protons, are shown. The spectrum of CRs at 
the given time and coordinates can significantly differ from the source 
spectrum. While the energy-independent diffusion leads to variation of 
the CR flux without change of the spectral form, the energy-dependent 
diffusion results in a modification of the particle spectrum. In Fig. 2 the 
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Figure 3. Gamma-ray emissivities in terms of E% x g 7 (_E 7 ) at different times t and 
distances R from a proton accelerator. The right-hand side axes shows the 7-ray 
fluxes, E? f x F 1 (E 1 ), which are expected from a cloud with parameter Mr,/d^ pc = 1. 
The thin and heavy curves correspond to times t — 1 3 and 10 5 yr, respectively. 
The fluxes at distances R = 10 and 30 pc are shown by solid and dashed lines. 
The power-law index, the total energy of protons, and the diffusion coefficient arc 
the same as in Fig. 2b. The curve shown by full dots corresponds to the 7-ray 
emissivity (and flux) calculated for local CR protons. In order to take into account 
the contribution of nuclei, all curves should be increased by a factor of r\ ~ 1.5. 
(From Aharonian and Atoyan, 1996). 

diffusion coefficient is assumed in a power-law form, D(E) oc E ^ above 
10 GeV, and constant below 10 GeV. The commonly used value for 
diffusion coefficient at 10 GeV, D 10 = D(10GeV), is about 10 28 cm 2 /s, 
however smaller values, e.g. D\$ = 10 26 cm 2 /s cannot be excluded, 
especially in the active star formation regions. 

The emissivity of the 7r°-decay 7-rays within 30 pc of relatively 
young, t < 10 yr, CR sources may significantly exceed the 7-ray emis- 
sivity introduced by the "sea" of GCRs, if the energy release in rela- 
tivistic particles per accelerator, like a supernova remnant or a pulsar, 
W50 > 0.01 (see Fig. 3). Under these circumstances the existence of 
massive gas targets like molecular clouds in these regions may re- 
sult in 7-ray fluxes detectable by EGRET, if (W 50 ■ M 5 )/dl pc ~ 0.1. 
In the case of energy-dependent propagation of CRs, large variety of 
7-ray spectra is expected, depending on the age of the accelerator, 
duration of injection, the diffusion coefficient, and the location of the 
cloud with respect to the accelerator. Therefore, at least part of the 
low-latitude unidentified EGRET sources, exhibiting a broad range of 
photon indices, T ~ 1.7 — 3.1 (Merck et al., 1996), could be explained 
as fortuitous location of CR accelerators in the proximity of massive 
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molecular clouds. A special interest could present also the mid-latitude 
EGRET sources which likely associate with the Belt Gould of massive 
stars and gas clouds (Grenier, 2000a; Gehrels et al., 2000). If this is the 
case, GLAST will be able to identify these EGRET hot spots, as well 
as to detect many more similar, but fainter (by one or two orders of 
magnitude) objects. 

The comparison of 7-ray fluxes from clouds located at different 
distances from an accelerator, may provide unique information about 
the CR diffusion coefficient D{E). Similar information may be obtained 
from a single 7-ray emitting cloud, but in different energy domains. For 
the energy-dependent propagation of CRs the probability for simulta- 
neous detection of a cloud in GeV and TeV 7-rays is not very high, 
because the maximum fluxes at these energies are reached at different 
epochs (see Fig. 3). The higher energy particles propagate faster and 
reach the cloud earlier, therefore the maximum of GeV 7-radiation ap- 
pears at the epoch when the maximum of the TeV 7-ray flux is already 
over. Meanwhile, in the case of energy-independent propagation (e.g. 
due to strong convection) the ratio F 7 (> 100MeV)/F 7 (> 100 GeV) is 
independent of time, therefore the clouds which are visible for EGRET 
at GeV energies would be detectable also at very high energies, provided 
that the spectral index of the accelerated protons a < 2.3. 



3. On the level of the "sea" of galactic cosmic rays 

The effective mixture of CRs from individual sources, during their 
propagation in the interstellar magnetic fields on time-scales r > 10 yr, 
determines the average density of CRs throughout the galactic disk - 
the level of the "sea" of GCRs. The expected small variation of the 
latter on large galactic scales does not, however, exclude non-negligible 
gradients of the CR flux on smaller (typically < 100 pc) scales, in 
particular in the proximity of young CR accelerators. 

It is generally believed that the local, i.e. directly measured at the 
Earth, flux of CRs (LCRs) gives a correct estimate for the level of the 
"sea" of GCRs. However, strictly speaking, this is an ad hoc assump- 
tion. In fact, it is not obvious that LCRs should be taken as undisputed 
representative of the whole galactic population of relativistic particles. 
In other words, we cannot exclude that the flux of LCRs could be 
dominated by a single or few local sources^, especially given the fact 

1 This statement is certainly true for the observed TeV electrons which suffer 
severe synchrotron and inverse Compton losses, and thus reach us, for any reason- 
able assumption about the diffusion coefficient, from regions no farther than a few 
hundred parsecs (Nishimura et al., 1980; Aharonian et al., 1995). 
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Figure 4- Gamma-ray emissivity in the Orion region. The experimental points are 
from Digel et al. (1999). The spectra of 7r°-decay 7-rays are shown by solid lines. 
The power-law (p-l) components at low energies, with the photon index T = 2.4, 
are shown by dotted lines. The dashed lines correspond to the total, " Tv°+p-l " 
spectra, (a) 7r -decay 7-ray emissivity calculated for the LCRs flux (curve 1), and 
for a 1.5 times lower flux (curve 2). The dashed curve corresponds to the sum of the 
latter and the low energy p-l component; (b) 7r°-dccay 7-ray emissivity (solid curve) 
calculated for the CR proton flux with power-low index index a = 2.1 and energy 
density w — 0.7eV/cm . The dashed curve corresponds to the superposition of the 
7r°-decay and p-l components. Note that the flux of the p-l component is assumed 
somewhat higher than in Fig. 4a. 



that the Solar system is located in a extraordinary region - inside 
active star formation complexes of the the so-called Gould Belt with 
essentially (more than factor of 3) enhanced core-collapse supernovae 
rate (Grenier, 2000b). Because of possible contamination of the "sea" 
of GCRs by nearby CR sources, we may expect significant deviations of 
both the spectrum and the energy density of LCRs from the spectrum 
and density of GCRs. Apparently, the "GCRs=LCRs" is an attractive 
hypothesis which however still lacks robust observational tests. Such an 
inspection can be uniquely done by observations of high energy 7-rays 
from GMCs. Although this method generally requires detailed spectro- 
scopic measurements for an ensemble of individual GMCs with known 
distances dk pc and masses M5, a reliable detection of 7-radiation from 
even a single under-luminous GMC, i.e. a cloud emitting 7-rays below 
the expected "standard" flux (dotted curve in Fig. 3), would imply that 
the CR flux we observe at the Earth is our local "fog" obscuring the 
genuine flux of GCRs. Below I argue that the 7-observations of the 
Orion complex contain evidence for such a dramatic conclusion with 
important astrophysical implications for the origin of CRs. 

The extensive study of high energy diffuse emission towards Orion 
by EGRET allowed an accurate measurement of the emissivity in the 
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Figure 5. Differential 7-ray emissivities in the Orion region. The curves 1, 2 and 
3 correspond to emissivities calculated for CR spectra with the same power-law 
index a = 2.1, but with exponential cutoff at 3 different energies, Eo — 100 GeV, 1 
TeV, and 10 TeV, respectively. The corresponding energy densities of CR protons 
are: w D = 0.55, 0.7, and 0.85 eV/cm 3 . At low energies an additional power-law 
component of radiation (cx E~ 2A ) is also assumed (dotted curve). The superpositions 
of 7r -decay and power-law components are shown by dashed curves. 



Orion region above 100 MeV, g 7 /4vr = (1.65 ± 0.11) x 10 -26 s _1 sr _1 , 
and, more importantly, a derivation of the differential emissivity in a 
broad, 30 MeV to 10 GeV energy region (Digel et al., 1999). These 
measurements, together with 7-ray emissivities calculated under dif- 
ferent assumptions about the CR flux and spectrum, are shown in 
Fig. 4. The emissivity of 7r°-decay 7-rays that corresponds to the local 
proton flux represented by Eq. (2), assuming an additional 50 per cent 
contribution from a-particles and other nuclei, is shown in Fig. 4a 
(curve 1). Although the predicted integral emissivity above 100 MeV, 
<7 7 /47r ~ 1.8 x 10~ 26 s -1 sr -1 , almost coincides with the measured emis- 
sivity, it is difficult to find a good fit to the derived 7-ray spectrum 
without assuming an additional low energy, e.g. a power-law component 
(dotted line). The latter with photon index T ~ 2.3 — 2.5 could be 
naturally attributed to the bremsstrahlung of primary or secondary (i.e. 
vr^-decay) electrons. It is seen from Fig. 4a that the superposition of 
the 'power-law' and '7r°-decay' contributions satisfactorily fits (dashed 
line) the derived experimental points, if we ignore the last, 1 to 10 GeV 
EGRET point. Note, however, that such a fit is achieved assuming a 
reduced, by a factor of 1.5, CR flux in Orion compared with the local 
CR flux. Perhaps, we should not overemphasize this difference, taking 
into account non-negligible systematic errors in the reported fluxes of 
local CRs, as well as uncertainties in the hydrogen column density used 
for derivation of the 7-ray emissivity. 



aharonian_gmc.tex; 1/02/2008; 16:24; p. 8 



Gamma Rays from Molecular Clouds 



9 



However, the discrepancy becomes really significant if we include in 
the fit the measured highest energy point at 1-10 GeV. This point ap- 
parently requires harder CR spectrum. The lack of 7-ray measurements 
above 10 GeV do not allow robust constraints on the power-low index, 
but, most probably, it should not exceed 2.3. In particular, Fig. 4b 
shows that a good fit to the EGRET data could be achieved assuming 
a flat spectrum of protons with power-law index a = 2.1, and energy 
density of particles below E < 1 TeV at the level of w = 0.7eV/cm 3 , 
which is by a factor of 1.7 less than the energy density of local CR 
protons given by Eq. (2). Note that the fit in Fig. 4b does not include 
the 7-ray flux produced by the "sea" of GCRs. This implies that the 
bulk of the observed 7-ray flux is contributed by particles accelerated 
within the Orion complex, with a significantly less (although unavoid- 
able) contribution from GCRs homogeneously distributed throughout 
the galactic disk, and freely passing through the molecular clouds. An 
extra contribution by GCRs around 1 GeV would lead to an over- 
production of 7-rays, unless we assume that the level of the "sea" 
of GCRs at energies ~ 1 — 100 GeV does not exceed, for the given 
7-ray flux uncertainties, approximately 1/3 of the directly measured 
flux of LCRs. Such a dramatic conclusion formally could be avoided 
speculating that CRs cannot freely penetrate into the dense clouds 
(see, however, Cesarsky and Volk, 1978). For example, assuming that 
the coefficient of reflectivity of CRs from a cloud decreases with energy, 
e.g. oc E~ K , with k ~ 0.6 — 0.7, one can fix the required hard spectrum 
of particles without invoking additional CR accelerators inside Orion. 

Both these possibilities generally do not agree with the conventional 
models of galactic cosmic rays. Nevertheless, a thorough analysis based 
on future comprehensive 7-ray data is obviously needed before any 
strong claim about a necessity for revision of the current basic concepts 
of origin and propagation of CRs. In particular, it should be noticed 
that the conclusion about hard CR spectrum in Orion is based merely 
on the reported 7-ray flux above 1 GeV. This indicates the importance 
of extension of accurate spectroscopic measurements of 7-radiation 
from Orion, as well from other individual 7-ray emitting clouds detected 
towards Cepheus, Monoceros, and Ophiuchus (Digel et al., 1999), to the 
energy region beyond 10 GeV. Also, it would be important to inves- 
tigate contributions of other 7-ray emission processes, in particular of 
the inverse Compton component contributed by relativistic electrons, 
albeit in dense gas regions the 7r°-decay radiation is believed to be the 
dominant 7-ray production channel. 

These measurements should allow an effective separation of contri- 
butions from the "sea" of galactic CRs and from the local CR acceler- 
ators. It is especially important to extend the measurements to higher 
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energies. If the "excess" flux at 1-10 GeV from Orion is indeed due 
to the local accelerator (s), we may expect continuation of this hard 
spectrum well beyond 10 GeV, depending on the high energy cutoff in 
the CR spectrum (see Fig. 5). On the other hand, the hypothesis of 
the energy-dependent reflection of galactic CRs from dense clouds (as- 
suming that GCRs=LCRs), predicts an essentially different spectrum 
of 7-rays - a flat part below a few GeV, with a quite steep power- 
law tail with T ~ 2.75, at higher energies. This part of the spectrum 
is produced by protons of the "sea" of GCRs which freely enter the 
cloud. Such standard high energy 7-ray tails should be observed from 
other local clouds as well. 

Large-scale CO surveys of molecular clouds by Dame et al. (1987) re- 
vealed two dozens of local GMCs within 1 kpc with M^/d^ > 10. The 
study of high energy 7-rays from these clouds with a broad-distribution 
of distances from ~ 100 — 200 pc (like the Taurus dark clouds or Aquila 
Rift) to 800 pc (like Cyg OB7 and Cyg Rift), which should be easily 
visible for GLAST, presents a great interest for derivation of spatial 
and spectral distributions of CRs in our local environment. Moreover, 
it is likely that GLAST will be able to extent the study of massive 
GMCs, and therefore to probe the flux of galactic CRs, beyond 1 kpc. 
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